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Gravitational lensing of the microwave background by the intervening dark matter mainly arises 
from large-angle fluctuations in the projected gravitational potential and hence offers a unique 
opportunity to study the physics of the dark sector at large scales. Studies with surveys that cover 
greater than a percent of the sky will require techniques that incorporate the curvature of the sky. We 
", lay the groundwork for these studies by deriving the full sky minimum variance quadratic estimators 

q ■ of the lensing potential from the CMB temperature and polarization fields. We also present a general 

technique for constructing these estimators, with harmonic space convolutions replaced by real space 
t products, that is appropriate for both the full sky limit and the flat sky approximation. This also 

extends previous treatments to include estimators involving the temperature-polarization cross- 
ed ■ correlation and should be useful for next generation experiments in which most of the additional 
information from polarization comes from this channel due to sensitivity limitations. 
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I. INTRODUCTION 



Weak gravitational lensing of the microwave background anisotropies offers a unique opportunity to study the dark 
matter and energy distribution at intermediate redshifts and large scales. In addition to producing modifications 
in the CMB temperature and polarization power spectra Q, 0, lensing of the CMB fields produces higher-order 
£^ ' correlations between the multipole moments 0, 3 • Quadratic combinations of the CMB fields can be used to form 
estimators of the projected gravitational potential, and therefore of the projected mass 0,0- The minimum variance 
quadratic estimator can in principle map the projected mass on large angular scales out to multipole moments of 
L ~ 10 2 and contains nearly all of the information in the higher moments of the lensed temperature field 0. 

Substantially more information lies in the lensed polarization fields allowing high signal-to-noise lensing reconstruction 



Of 

O ' and extending the angular resolution out to L ~ 10 3 (lCf. 
+3 . Lensing reconstruction techniques involving the polarization fields have previously only been developed for small 
surveys where the sky can be taken to be approximately flat. Since lensing is intrinsically most sensitive to the 
projected potential at L < 10 2 or several degrees on the sky, a treatment incorporating the curvature of the sky 
is desirable. In fact it is necessary for its application in removing the lensing contaminant to gravitational wave 
| polarization 0, 0, across large regions of the sky. 

We present a concise treatment of the effect of gravitational lensing on CMB temperature and polarization harmonics 
in Sect.[n] We construct the full sky quadratic estimators of the lensing potential and compare their noise properties 
to that for the flat sky expressions in Sect. IIIII We provide an efficient algorithm for the construction of all estimators 
in Sect. llVl We summarize some useful properties of spin- weighted functions in Appendix 1X1 Finally, we derive the 
flat sky limits of the estimators and draw the connection to results in [lCj in Sect. iBl 



II. CMB LENSING IN MULTIPOLE SPACE 



In this section, we give a ped agogica l but concise derivation of the lensing effect on the CMB temperature and 
polarization fields on the sphere |L3L Il4| . We emphasize the connections between the formalism using spin- weighted 
spherical harmonics |15| and a tensorial approach [Tfjl ] which will be useful for the lensing reconstruction in the 
following sections. 
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The temperature perturbation is characterized by a scalar function 8(n) = AT(n)/T, whose harmonic transform 
is given by 

e(n) = 5>P5T(n)- (1) 

im 

The polarization anisotropy of the microwave background is characterized by a traceless, symmetric rank 2 tensor, 
which can be represented as (e.g. [T^ l 

Vij = + 2 A(n)m,mj + -2A(h)m l m : j 1 (2) 
where we have defined the complex Stokes parameters ±2A according to 

±2 A{h) = Q{h)±iU{n). (3) 
The spin projection vectors are given with respect to the measurement basis (ex, e 2 ) by 

m = -j= [ex + ie 2 ] , (4) 

m = [ex - ie 2 ] , (5) 

and form an eigenbasis under local rotations of basis vectors (see Appendix |A"|) . In spherical polar coordinates, ex = eg 
and e.2 = e v . Under a local, right-handed rotation of the basis (ex, e 2 ) by an angle ip, the complex Stokes parameters 
±2A(h) acquire a phase e T21 ^ . They act as spin-2 functions, with a corresponding harmonic transform in terms of 
spin-weighted spherical harmonics |17| given by |l5j 

±2 A(n) = £ ±2 ^ ± ^™(n). (6) 
i m 

A lens with a projected potential 0(n) maps the temperature and polarization anisotropies according to |3.l3Lll8| 

8(n) = 6(n + V^(n)) 

= e(n)+V^(n)V 4 9(n)+O(0 2 ), (7) 

Vij(n) = T> lj {n + X7cf>(n)) 

= T t] {n) + V k 4>{n)V k T l3 (n)+0{4> 2 ), (8) 

where tildes denote the unlensed fields. In the case of a weak gravitational field under consideration, lensing potential 
cj> is obtained by a line-of-sight projection of the gravitational potential, 

where 7/ is the conformal time, rj s is the epoch of last scattering and x is the angular diameter distance in comoving 
coordinates. 

Taking the harmonic transform of Eqn. J2J, one readily shows that the the change to the temperature moments 
5&l n ee Of 1 - 6J™ are given by 

SQ r * EE^©™' 7 ^ m '> ( 10 ) 



with IJ^fl 1 " 1 denoting the integral 



LM I'm' 



Iltv m ' = / dnY^ViY^VYf ' . (11) 



The integral can be performed analytically using the relation |2 

3 



dn Sl Y^(n) S2 Y^(n) S3 Y^(n) 



1/2 



-in 



H 2h + 1 



h h h \ ( h h h 
-sx -s 2 -S3 ) I mi m 2 m 3 



(12) 
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to yield 



jraMra 



{-iy 



l L V 
—m M m' 



oFi 



ILV 



with the definition 



Fi Ll > = [L(L + 1) + + 1) - 1(1 + 1)] 



(2L + l)(2i + l)(2Z' + l) (ILV 
16tt Us Ts 



(13) 



(14) 



The multipole expansion for the polarization fields proceeds by noting that ±iA(n) are the spin ±2 components 
of the polarization tensor "Py. Since the contraction with the spin projection vectors m l m? projects out the spin 2 
piece of a symmetric tensor, the change in the complex Stokes parameters is given by 



5[ +2 A{n)} = m l m j SVi 



VkfiAh) [V fe 0(n)]. 



(15) 



The expression for the contribution to _2^4(n) is obtained by replacing m % by rri 1 in the above. We denote the 
product m'm^V kPij using a spin-gradient derivative Di (see Eg. IA12| . and write the lensing contribution as 



S[±2A(h)} «£>V(n)A[± 2 i(n)]. 



(16) 



This relationship was given in [l3] | with the shorthand convention Di — > correspondin g to the action of covariant 
derivatives on the spin components of symmetric trace free tensors given in Eqn. ijAllfl |l4l Il7| . Expanding <p and 
±2^4 in spin- weighted spherical harmonics and evaluating the inner product of their gradients using Eqn. ijATSf . we 
obtain the lensing corrections 



0[±2A l \ « Z^ ( PL ±2 A l ±2hU' > 



LM I'm' 



where we define 



rraMra' 



±2 1 ILV 



(-iy 



l L I' 
-m M m' 



±iFiu 



We will be interested in the lensing expressions for the rotationally invariant combinations 



Bl ~ 27 



+2A\ 



-iA 



2^1 J , 



(17) 

(18) 

(19) 
(20) 



which are the curl- free ("E-mode") and gradient-free ("B-mode") components of the polarization field. ^From the 
expressions (|10|) and (|17J) . we find the general expression for a lensed multipole moment to be 



LM V 



I I' 



L 



m —m! —M 



(21) 



where X™ may be multipole moments of O, E, or B, and 



l + (-l) L +'+ r 



ft, 



IV IL — 



l-(-l) 



L+l+l' 



2? 



(22) 



ensure that the associated terms are nonzero only when L + I + I' is even or odd, respectively. X denotes the parity 
complement of X, i.e. 6 = 0, E — —B, B = E. 
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III. QUADRATIC ESTIMATORS 

Lensing of the CMB fields mixes different multipoles through the convolution 121fl , and therefore correlates modes 
across a band determined by the power in the deflection angles [13j . The unlensed CMB multipoles AT™ are assumed 
to be Gaussian and statistically isotropic, so that the statistical properties are characterized by diagonal covariances 
or power spectra 



(23) 



The assumption of parity invariance implies that Cf B = Cf B = 0. The lensing potential is also assumed to be 
statistically isotropic so that 



l(l + l)(4>r<l>™) = S„,5 mm ,Ci 



(24) 



where we have multiplied through by 1(1 + 1) to reflect the weighting of deflection angles. 
It follows then that the lensed multipoles are also statistically isotropic with power spectra 



(X ; m *X'™ ^ _ Su'6 mm >Ci 



XX' 



(25) 



Since XJ 11 denotes the measured multipoles, the power spectra contain all sources to the variance, including detector 
noise. Detector noise will be taken to be homogeneous, with power spectra given by [2l| 



j/81112 



a 



EE i 



■ = c, 

noise 1 



BB I 



(j^y 

V Tqmb ) 



j/81ii2 



(26) 



where Ae and Ap characterize detector noise, and 6*fwhm is the FWHM of the beam. We employ the specifications 
of a nearly ideal reference experiment, with Ae = l^K-arcmin, A P = %/2>K -arcmin, and 0fwhm = 4' (see 0] for 
an exploration of noise properties). 

If we instead consider an ensemble of CMB fields lensed by a fixed deflection field, the multipole covariance acquires 
off-diagonal terms and becomes 



/ m t rn 



lens 



M 



LM 



I I' 



L 

-M 



ft 



a ±M 
ILVfh i 



(27) 



where the subscript on the average indicates that we consider a fixed lensing field. f{* Ll , are weights for the different 
quadratic pairs denoted by a, given by 



fhLl 2 = s a Fi lL l 2 £hi 2 LCi 2 + f3i 1 i 2L C l 



y^iab a s~iat 



ah 



(28) 



where s a and Sb are the spins of the a and b fields respectively. Specific forms for the six quadratic pairs are given in 
Table H| 

/a 
hLl 2 



ee 

@E 
EE 
QB 
EB 
BB 



Oh oFi 2 Lh + Cf 2 ~ oFi lL i 2 
Cf^iFi 2th + Cf 2 E F hLh , even 
Cf 1 E 2 Fi 2Lll + Cf 2 E 2 F lL , 2 , even 



iC, 
iF, 



t L E 2fl 2 Ll 1 , 



2 r l 2 Lh 



c. 



odd 

2 F, 



l 2 2-t i l 1 Ll 2 , 



odd 



riBB r . fiBB r 

Cjj 2fi 2 Li 1 +Ci 2 2l<i 1 Li 2 , even 



TABLE I: Functional forms for f t a lLh 
or odd, respectively. 



"Even" and "odd" indicate that the functions are non-zero only when L + I1 + I2 is even 



Because has a zero mean, the off-diagonal terms of the two-point correlations (a™ 6™ ) taken over a statistical 
ensemble would vanish. However, in a given realization, we can construct an estimator for the deflections as a weighted 
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sum over multipole pairs, and find weights that minimize the variance of the estimator. We write a general weighted 
sum of multipole pairs as 



A°, 



1) , , V mi m 2 —M I 



(29) 



where a™ and 6™ are the observed CMB multipoles, a denotes the specific choice of a and b, and the sum includes 
the diagonal (h — I2, mi = —m 2 ) pieces. 

^From expression Ij27(l for the average over a fixed lens realization, 



■a M \ 



Al 



lens y/L(L+ 1) 

where we have used the relations 



2L 



hh 



E 



1 



h h L }( h h V \ 

mi m 2 M I \ mi m 2 M' ) 2L + 1 



Sll'Smm', 



771 \ 



I I L 
m —m 



21 + 1 
2L + 1 



'Lo- 



rn 

(31) 
(32) 



The diagonal terms in Eq. I|29|l only contribute to the unobservable monopole piece and we hereafter implicitly 
consider L > only. The normalization is set by the condition 



M 



(33) 



to be 



A a L = L(L + l)(2L + 1) E 9? lh (L)f? iLl2 



(34) 



.hh 



We derive the minimum variance estimator by minimizing the Gaussian variance {d a f *d a l I ) with respect to gf lh {L) 
and find that 



9hh(L) 



Note that for a = b, 



h h J l\Lh V / h h JhLh 

(Jaa (Jaa (jbb (jbb (C a ^C a ^ )^ 

h h h h ^ h h ' 

£OL* 
JhLh 



and for Cf = (e.g., for QB or EB), 



9hi 2 (L) 



9hh{L) 



2Cf a Cf a ' 



■ T h 



u 2 



rjaa^jbb ' 
h 1-2 



The Gaussian noise covariancc 



is given by 



N 



a/3 



L{L + 1){2L- 



T) E {dhW [qrctf^w + (-i) L+ ' i+ ' 2 Qfcf 2 c 3f 2il (^)] } 



(35) 



(36) 



(37) 



(38) 



(39) 



hh 



with a = (ab), (3 = (cd). For a — f3, the above reduces simply to N£ a = Al- 
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FIG. 1: Deflection and noise power spectra for the quadratic and minimum variance estimators, assuming the noise properties 
of the reference experiment (Ae = 1/iK-arcmin; Ap = -v/S^fv-arcmin; fl FWHM = 4') and a fiducial ACDM cosmology with with 
parameters Q, c — 0.3, fif, = 0.05, Qa = 0.65, h = 0.65, n = 1, 8h = 4.2 x 10 -5 and no gravitational waves. 



Following the treatment of the flat sky case in |lOj |. we combine the measured quadratic estimators to further 
improve the signal to noise by a forming minimum variance estimator 



<T v f = ^w a (I)d a f, (40) 



with weights and variance given by 



w a (L) = JVr£( N ^) a/J ' (41) 



N l V = 7^7^- ( 42 ) 

(43) 

We will hereafter ignore contributions from the BB estimator, since the primordial contributions to the i?-mode power 
spectrum is expected to be small on scales where the lensed multipoles are employed. 

We plot the noise power spectra for the five estimators, as well as the minimum variance estimator, in Fig. ^ 
assuming the noise properties of the reference experiment. 



IV. EFFICIENT ESTIMATORS 



The quadratic estimators involve both filtering and convolution in harmonic space. It is useful in practice to 
express the convolution as a product of the fields in angular space. The estimators can then be constructed using fast 
harmonic transform algorithms p3 , |23| . To simplify the construction of the estimators we will assume Cf B « Cf E 
as is appropriate for the standard cosmology. Aside from the EB estimator, derived in [lfj under the flat sky 
approximation, the angular space estimators involving polarization are new to this work. 

Generalizing the construction in @ for the 00 estimator, consider the fact that lensing correlates the (lensed) 
temperature and polarization fields to the their (unlensed) angular gradients. We show in Appendix ^ that the 
all-sky analog to the gradient operation on a spin-s field is di — > Z?,;. The quadratic estimator is then built out of the 
general operation on two fields X(h) and Y(n) 



P[X(n),Y(n)} = -D l [X(n)DiY(h)} 



(44) 
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The properly normalized estimators then take the form 



s/L{L + 1) 



dny i M, (n)e a (n). 



(45) 



where 



e ee (n) 


= P 


A e , o^ee 


] ; 




e eE (n) 


1 

~~ 2 


(P[ +2 A £ ,_ 


2Ae_E] - 


hP[cc}) + P[ A e , A E e 


e eB (n) 


1 

~ 2 


(P[ +2 A iB , 


-2Aqe\ 


+ P[cc]) , 


e EE (n) 


1 

~~ 2 


(P[ +2 A £ ,_ 


2Aee] - 


hP[cc]) , 


e EB (h) 


1 

~~ 2 


(P[ +2 A. iB , 


-2Aee] 


+ P[cc]) , 



(46) 

where cc denotes the operation with the complex conjugates of the fields and the filtered fields themselves are given 
by the general prescription 



±sAx (n) = Y /? ixXr±sYr(n), 



± s A X Y(h) 



c<xy 

E^i 
(J xx 



(47) 



We omit the a — BB estimator under the assumption that the unlensed P-power is small at high multipolcs 

L ~ a L 



It is straightforward to verify that all of the estimators are the same as the harmonic space ones d al £ — d a ^ with 



At = A L except for QE. Here the weights on the multipole combination are 



n 0E 



fir) 
■111 



eE 

Ll 2 



r<9&(jEE 
h h 



(48) 



and are slightly non-optimal compared with the minimum variance weighting. Furthermore Nl =/= Al and they 
must be calculated separately. However a direct calculation of the noise spectrum through Eqn. 139|) shows that the 
differences are less than 1%, and essentially indistinguishable from the minimum variance estimator (see Fig- - 

These estimators may therefore be used in place of a direct multipole summation for efficient lens reconstruction. 
The gradient operations in Eqn. (|44() are efficiently evaluated in harmonic space since their action on spin harmonics 
simply raises and lowers the spin index in accordance with Eqn. (|Af 71) . 



V. DISCUSSION 



Counterintuitively, the gravitational lensing of the CMB temperature and polarization fields is a small scale mani- 
festation of the very large scale properties of the intervening mass distribution. It therefore requires very challenging, 
high angular resolution (< 10') but wide-field surveys (> few degrees) to exploit. We have provided expressions for 
quadratic estimators of the lensing potential valid on the entire sky, as well as the expected noise covariances for 
the estimators. As expected, on small angular scales (L > 100), the flat sky approximations differ from the full sky 
expressions by less than ~ 1%, indicating that the flat sky approximations is adequate. This regime is however not 
where the signal-to-noise peaks. 

We have also provided a practical means of implementing these estimators using fast harmonic transforms, either 
with spherical harmonics or Fourier harmonics, to perform the required harmonic convolutions and filtering. We have 
shown that even the approximate QE estimator has a noise performance that is essentially indistinguishable from the 
minimum variance estimator. These techniques should provide a means to study the impact of real world issues such 
as finite-field, inhomogeneous noise, and foregrounds on the science of CMB lensing. 
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FIG. 2: Fractional difference between the approximate and the minimum variance QE estimators, defined as 8Nl/Nl = 

/^•approx. _ 7V"mv)/jV™ v . 

APPENDIX A: TENSOR REPRESENTATION OF SPIN WEIGHTED FUNCTIONS 

We clarify the relation between spin-s functions and tensor quantities on the sphere, and derive the relation between 
spin raising and lowering operators and covariant derivatives on the sphere. 

Suppose we construct an orthonormal basis (ei(n), e.2(n)) at each point on the sphere, with n denoting the outward- 
facing normal vector. We define a local rotation as a right-handed rotation of the basis vectors (ej., e 2 ) by an angle ip 
around the vector n, so that the new basis vectors [ey, ey) are related to the original vectors by the transformation 



e.y = cos ipei + sm.rpe.2i 
&2' = — sin ip e-i + cos ip $2- 



(Al) 



A function s f(n) is said to carry a spin-weight s if, under the rotation (|A1|) . the function transforms as s /(ti) 
e~ ls ^ s f{n). This convention conforms to ^5|, and defines rotations in a sense opposite to that in [TtI |24|. 



We define vectors m and m with respect to the basis (ex, e 2 ) according to 

1 



which have the property that 



m = 



m = 



m ■ m 
m ■ fn 



V2 
1 

71 



ex + ie 2 \ 



ex - ie 2 \ 



m ■ rn 
1. 



(A2) 
(A3) 

(A4) 



Given a vector field v(n), it can easily be shown that the quantities v ■ m and v ■ fh transform as spin 1 and — 1 
objects, respectively, so that m and m act as spin projection vectors. More generally, given a rank-s tensor T^...^ , the 
quantity T ix is m 11 ■ ■ ■ m> transforms as a spin-s object, since under the rotation (|A1|) . each factor of m %n contributes 
a phase e~ ls ^ . 

The spin-s functions therefore also provide a complete basis for the totally symmetric trace-free portion of a 

rank-s tensor 



3 /(n)m i:1 • • • m ls + _ s /(n)m il ■ • ■ m l 



(A5) 



where the trace-free condition refers to the vanishing under contraction of any two indices in the tensor. For example, 
the polarization tensor can be written as 



(A6) 
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Covariant differentiation of such a tensor is related to the raising and lowering of the spin weight : 

VfeT^...^ = [d fcs /(n)] ■ rh lg + s /(n)V (fc m il ■••my 

+ [d k - s f(n)}m il ---rn ls +- s f(ti)X7 {k m. ll ---m is) . (A7) 

We evaluate the covariant derivatives Vifhj etc. explicitly in the spherical basis with coordinates (0,(p), yielding 

Vgrhg = Wgfh v = 0, 



\7 v mg = —=cos0, 
V2 

V^m^ = -^sin#cos#, 

with those for m given as complex conjugates of the above. Using these, it can be shown that 

0, 



m- 7 V 'itrij 



m J V irrij = 
—rh 3 V 'itrij 



Ji , 



where 



Jg = 0, J v — i cos ( 



in spherical coordinates. The covariant derivative of T^...^ is therefore given by 

VfeT^...^ = [D k s f(n)] rh^ ■ ■ ■ rh is + [D k _ s /(n)] m h 
where we define the spin-dependent gradient operator as 

Di = di + sJi. 



■mi 



(A8) 

(A9) 
(A10) 
(All) 
(A12) 



A covariant derivative operating on the spin-s piece of a tensor is equivalent to a gradient operation Di on its 
spin-s weighted representation. As an example, the components of the covariant derivative of the polarization tensor 

m/m'YiT,, = D k [ 2 A{h)], 

m'm/YiT,, = D k {- 2 A(n)}. (A13) 

The gradient operator D k is related to spin raising and lowering operators. Using the expressions IjAlOfl and 
expressing the operator D k in the (m, ffi) basis, we obtain the desired relations 

1 



A [»/(«)] = ~-j= {[&/(«)] m i + dsf(n) m 4 j. 



(A14) 



By virtue of the rotational properties of (m, m), the ladder operators d and d, defined by ^3,13 

6.f(6,<p) = -sin s # — +icsc6— sin- s 9 s f(9,<p), (A15) 

8 /(6> )¥ >) = -shU^ zcsc^— sm s 9 s f(9,<p), (A16) 

raise and lower the spin weight by 1. For example, the gradient operation on the spin-s spherical harmonic yields 



DiUvn = ~ (hi - s)(i + s + i^s+iv 



m, 



[(l + sW-S + l)] 1 ' 2 .-!^**) 



(A17) 



Note that the inner product of two gradients 

[D* ai fi(fl)][Dis a Mn)] = \ {[$.Ji(n)] [^ a /a(n)] + [6 s Ji(n)] [3 S2 .f 2 (n)}} (A18) 

leaves the total spin-weight of the product unchanged. 

Inverting the relation (|A14(I . we obtain the ladder operators in the tensor representation 

d„f(n) = \ Im rn ■m i V,/, , . (A19) 

S a f(h) = \ 2m m' m' V,/, , . (A20) 

for s > 0, and with m 4 ™ replaced by fh ln for s < 0. This relationship was first proven in ^tJi albeit with a different 
sign convention. 
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APPENDIX B: FLAT-SKY APPROXIMATION 



The all-sky estimators derived in Sect. IIIII reduce to the flat sky estimators, based on Fourier harmonics of the 
fields in the small angle limit. Here we explicitly show this correspondence. 

The full sky harmonics in multipole space (I, m) are related to the flat sky harmonics in Fourier space I = I cos ipix + 

zsin^yby mm 



X(l) 

Rewriting Eq. (|29(l using the above, 



4tt 



21 + 



- i~ m X, 



m imipi 
i e 



M 



i) e rL > > i 
(2/i + 1)(2Z 2 + 1) f dLp x dip 2 



hmi l 2 rri2 



h h L 
n%i m 2 —M 



y 4tt(2L + 1) 7 2tt 2tt 
To go further, we can utilize the approximation [T: 



e -* mi¥ ' 1 e- !im2V2 a(«i)6(i2)- 



(Bl) 



(B2) 



h L h \ M x L Z 2 \ I cos2(i^ 2 - ^J, L + Zi + Z 2 = even, 
2 O-2/lo 0/ I isin2(p j2 - ^J, L + h + l 2 = odd, 



(B3) 



with the trigonometric functions defined through the cosine and sine rules, and the 3-j symbol on the rhs for the odd 
case represents a continuation of the analytic expression for the even case 



la lb lc 





= (-1)*/ 2 



(1)1 



ti-l*Mk-h)l(h-ie)l 



(i-2i a y.(i-2i b y.(i-2i c y. 



(i + iy. 



1/2 



(B4) 



where I — l a + lb + l c - 

Furthermore, in the limit 1%, Z 2 , L 3> 1, 



1 



[L(L + 1) + l 2 (l 2 + 1) - h(h + 1)] « L ■ Z 2 , 



(B5) 




FIG. 3: Fractional differences 8Nl/Nl between the noise in the flat sky and full sky estimators, calculated for the reference 
experiment. 
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and so we may absorb the geometric factors in a FliH? a $ 



with 



. / (2L+l)(2Z 1 + l)(2/ 2 + l) hl 2 L\ p 
sF hL i 2 ~\l ^ I \ s F hLh , (B6) 



1, s = 0, 

= L ■ l 2 x ^ cos2(^ 2 — </5i 1 ), s = 2, £ + + Z 2 = even, 
i sin2(i^; 2 — (fi t ), s = 2, L + li + l 2 = odd. 



(B7) 



The approximations for f^ L i, and gf^ ^L) can likewise be written as 



TOi ^ / V"— i "A""t ' ^yV"-"^ 1 *y I "1 "A I /*CK 

/jiifa - \/ ^ ^ Jt h Ll a , 

n a m „ / (2£+l)(2Z 1 + l)(2Z 2 + l) ( l x h L\ rn , , 

ffJibW ~ V 4^ I o o JdhiAL), (B8) 

where /;" L ; 2 and gf^{L) are defined as the unbarred quantities with sF^l^ replaced by s Fi 1 [ J i 2 . We will also utilize 
the relation between plane waves and spherical harmonics, 

e ihn, _ ^ i^yrn. e -^Vi . (B9) 

m 

Using Eq. l|B8Jl to rewrite Eq. <jB2fl . and applying relations (|12|) and l|B9|) . the estimator becomes 



/L(TTT)(2TTi)^ 



4tt 



Taking li,l 2 ,L^> 1, the above reduces to 

d«(X) « ^ J J0^g? ih ( L )a(h)b(h), l 2 = L -h (Bll) 

with g" il2 (L) corresponding to the filters F a (li, l 2 ) in [Tc}- The normalization reduces to the flat sky expression in 
|10| in a similar fashion, by using the approximations l|B8|) to relate the full sky quantities to trigonometric functions 
on the flat sky. 

It is simple to show that the efficient all-sky estimator in Eqn. (|45l) reduces to efficient flat sky estimators with 
the replacements Di — ► di in Eqn. I|44|) and the spherical harmonic transform in Eqn. (|45ll with a Fourier transform. 
Under the assumption that Cf B <§C Cf E , they again reproduce the properties of the minimum variance quadratic 
estimators in Eqn. (|Bllf) and allow fast Fourier transform techniques to be employed in their construction. 

Fig. shows fractional differences between the noise in flat sky estimators derived in ^(| and the noise in full 
sky estimators, defined as 5N£/N£ = (7V^ fiat ^ — j /v^ ful ^)/7V^ ( ' full ' ) . Because most of the information comes from 
multipole pairs at high multipole moments, the flat sky expressions deviate at less than ~ 1% for L > 200, mainly in 
the direction of overestimating the noise. 
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